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Abstract 

A simple set of algebraic equations is derived for the exact low-temperature thermodynamics of 
one-dimensional multi-component strongly attractive fermionic atoms with enlarged SU(N) spin 
symmetry and Zeeman splitting. Universal multi-component Tomonaga-Luttinger liquid (TLL) 
phases are thus determined. For linear Zeeman splitting, the physics of the gapless phase at low 
temperatures belongs to the universality class of a two-component asymmetric TLL corresponding 
to spin-neutral iV-atom composites and spin-(iV — l)/2 single atoms. The equation of states is also 
obtained to open up the study of multi-component TLL phases in ID systems of iV-component 
Fermi gases with population imbalance. 

PACS numbers: 03.75. Ss, 03.75.Hh, 02.30.1k, 05.30.Fk 
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A significant feature of one- dimensional (ID) many-body systems is the universal low 
energy physics of the gapless phase described by a Tomonaga-Luttinger liquid (TLL) [lj], 
as recently revealed in experimental measurements on the thermodynamics of spin ladder 
materials |2|. In such systems magnetic fields drive phase transitions between gapped and 
gapless phases. In the gapless phase spin excitations (spinons) carrying spin- 1/2 give rise to 
universal TLL thermodynamics. 

On the other hand, exquisite control over the effective spin-spin interaction between cold 
atoms has provided a new opportunity to rigorously test spin liquid behaviour via trapped 
fermionic atoms with higher spin symmetry 0, Q, E]. In particular, fermionic alkaline- 
earth atoms display an exact SU(N) spin symmetry with N = 21 +1 where / is the nuclear 
spin [8]. For example, I = 9/2 for 87 Sr and I = 5/2 for 171 Yb. Such fermionic systems 
with enlarged SU(N) spin symmetry are expected to display a remarkable diversity of new 
quantum phases due to the existence of multiple charge bound states. 

We derive the universal thermodynamics of the gapless phase in strongly attractive ID 
fermions with iV-component hyperfine states. We thus find that population imbalances con- 
trolled by Zeeman splitting can be used to explore and control multi-component TLL physics 
in ID interacting fermionic gases with SU(N) spin symmetry. The universal crossover from 
the regime of a relativistic multi-component Luttinger liquid to a nonrelativistic quantum 
critical regime is determined from the specific heat phase diagrams. 

In principle, the thermodynamic Bethe ansatz (TBA) can provide the exact thermody- 
namics of such systems. However, the TBA involves an infinite number of coupled nonlinear 
integral equations which hinders access to the thermodynamics from both the numerical and 
analytical points of view j^], [h]] . To overcome this obstacle, a method was proposed to ana- 
lytically obtain thermodynamic quantities of the two-component attractive Fermi gas jll ]. 
We build on this approach to reduce the multiple charge bound state TBA for attractive 
ID fermions with arbitrary SU(N) spin symmetry to a set of simple algebraic equations to 
provide universal and exact low-temperature thermodynamics of the model. 

The model. We consider a system of Nf interacting fermions of equal mass m which 
may occupy iV possible hyperfine levels i = 1, . . . , N) labeled by iV isospin states and 
constrained by periodic boundary conditions to a line of length L. The Hamiltonian [12I, [3] 
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IS 

h 2 Nf d 2 

H = -^Y^^2 + 9m Yl S(xi- Xj )+E z (1) 

i=l 1 l<i<j<Nf 

with Zeeman energy E z = YliLi ^Wz^b-i Here Nj is the number of fermions in state 
\i) with Zeeman energy e z (or say, a chemical potential //j) determined by the magnetic 
moments [i B and the magnetic field B. The spin-independent contact interaction g^ remains 
between fermions with different hyperfine states so that the spins are conserved, i.e., Nj 
with i = 1,...,N are good quantum numbers. The coupling constant gm = h 2 c/m with 
interaction strength c = — 2/aio determined by the effective ID scattering length am. For 
simplicity, we choose the dimensionless units of h = 2m = 1 for numerical calculation 
and use the dimensionless coupling constant 7 = c/n with linear density n = Nf/L. The 
Hamiltonian ([1]) exhibits U(l) x SU(N) symmetry where U(l) is associated with the charge 
degree of freedom and SU(N) with the N hyperfine spin states. 

Ground states. The energy eigenspectrum is given in terms of the fermion quasimomenta 
[hi] satisfying the Bethe ansatz equations (BAE) [12I, Q by E = Yuf=i kj- For attractive 
interaction, the BAE allow charge bound states and spin strings. The SU(N) symmetry 
acquires N — 1 kinds of charge bound states which can be viewed as composites of r-atoms 
with total spin s, where s = r(N — r)/2 and r = 2, . . . , N. In order to simplify calculations 
in the study of magnetism, we rewrite the Zeeman energy as Ez = — Yli^ 1 HiNi where 
the independent parameters Hi with % = 1, . . . , N — 1 characterize the chemical potentials 
for A r i unpaired fermions and iVj charge bound states of z-atoms. Here we set up explicit 
relations among the H^s and arbitrary (nonlinear) Zeeman splittings Aj + xj = e^ 1 — e l z 
via the relations A i+li = —H^i + 2Hi — H i+1 with H N+ i = 0. We shall show that equally 
spaced (linear) Zeeman splitting, i.e., A a+la = H for a = 1, ... , N — l, drives a smooth phase 
transition from bound states of iV-atoms into a normal Fermi liquid at zero temperature, see 
Fig. dtA). Nonlinear Zeeman splittings may trigger spin-s charge bound states as illustrated 
in Fig. [U(B), where the magnetic moments of paired states is s = N — 2 [4J]. 

For strong attraction (I7I ^> 1) these charge bound states are stable and the system 
is strongly correlated. The corresponding binding energies of the charge bound states are 
given by e r = h 2 c 2 r(r 2 — l)/(24m). Solving the BAE with strongly attractive interaction, 
the ground state energy E£° per unit length in the thermodynamic limit is given explicitly 
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by (in units of h 2 /2m) 



N .2^3 / - x N 



3r V c c 2 , 

r=l ' r=2 



Here n r = N r /L and we have denned the functions 



r-l N , „/ _s N 



4n i 0(r - 2) ^E 4r^(iV -r-l 

where 6*(x) is the step function. 



^ ^ r(i + r - 2j) ^ r(i - r) 

j=l i=j V J7 i=r+l v 2 



The general result ([2]) for arbitrary iV unifies and extends known results mm 1 161 ] for 



isospin 5 = 1/2, 1 and 3/2. Furthermore, we find that the low-lying excitations are described 
by the linear dispersion relations e r (k) = v r {k — k F ), where v r ~ ^^(1 + iff A r + ^A 2 ) with 
r = 1, . . . , N are the velocities for unpaired fermions and charge bound states of r-atoms. We 



denote the corresponding 



to the universal form 



13, 



ermi momentum by k r F . These dispersion relations lead naturally 



181 ] for finite-size corrections to the ground state energy 



hC N 

E(L,N f )-LE^^-^J2v r . (4) 

r=l 

The central charge C — 1 for U(l) symmetry. Here the universal finite-size corrections (T4|) 
indicate the TLL signature of ID many-body physics [l|. 

Quantum phase transitions. At finite temperatures, the equilibrium states are determined 
by the minimization of the Gibbs free energy {20!], which gives rise to a set of coupled 
nonlinear integral equations - the TBA equations 2l|. The TBA equations provide not 
only finite temperature thermodynamics but also the band fillings with respect to Zeeman 
splittings and chemical potentials which may be conveniently used to analyze quantum 
phase transitions at zero temperature. From the TBA equations, we find that complete 
phase diagrams at zero temperature are determined by the independent external field-energy 
transfer relations 

H r = ±r(N 2 -r 2 )c 2 + r(^-^ N ), (5) 

for r = 1,...,N — 1. The effective chemical potentials /i^ = \i + — + ^ T Z 1 ^ ^ with 
r = 1, . . . , N are given by 



M n 2 ( 2A r 3A 2 \ B r -I 3B r ■ A 



7T 2 r 2 V Id c 2 |c| c 2 



3 1 + ^ + ^ + ^ + (6) 
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characterizing Fermi surfaces of stable charge bound states and unpaired fermions. Here 
A = {At, . . . , A N ) and B r = (Bl,..., B?) with 



B. 



3r£ 



£ 

.i=i 



B(r-j) 



(r + £-2j) 



+ 



6(r-e- 1) 



r — 



(7) 



for r, £ = 1, . . . , N. I is a unit vector. The energy transfer relations (jSJ) can be used to 
show that linear Zeeman splitting can only lift SU(N) symmetry to U(1) N symmetry, recall 

Fig.m 

For pure Zeeman splitting, if the external field H is less than a lower critical field H&, 
a molecular superfluid phase of spin-neutral bound states forms the ground state. We find 
the value 

(iV+l)n 2 7 2 27r 2 n 2 f, 16Q 
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N*(N - 1) 

2og 2 



i + 



+ 



3N 2 \j\ 
16Q 



(8) 



3JV(iV-l)|7| 7 2 A^ 4 -f 2 N 3 (N-l)_ 

at which the spin gap is diminished. Here Q = J2f=i '■ ® n ^ ne other hand, when 
H > H C 2, where 



2^,2 



H, 



c2 



(N + l)n 2 i 2 ^ 2tt-ii 



1 - 



3N(N- 1)| 7 | 



(9) 



6 ' (JV-1 

the system is fully-polarized into a normal Fermi liquid. For the intermediate regime 
H c i < H < H C 2 spin- neutral bound states of iV-atoms and unpaired fermions coexist. 
The magnetization gradually increases from m z = to its normalized value m z = 1 as the 
field increases from H c % to (see the solid line Fig. [2]). More subtle phases may be ex- 
plored using the energy-transfer relations (jSJ) by controlling the nonlinear Zeeman splitting 
parameters, e.g., the BCS pairing phase (see Fig. [TJ by setting A 2 i = hi and A a+ i a = h 2 
for a = 2,N — 1. All phase transitions are of second order with a linear field-dependent 
magnetization in the vicinity of critical points due to the condition of fixed total particle 
number. 

Universal thermodynamics. At high temperatures T ~ e r (setting the Boltzmann constant 
to unity), thermal fluctuations can break the charge bound states while spin fluctuations 
lead to an effective ferromagnetic spin-spin interaction coupled to each Fermi sea of spin-s 
charge bound states and unpaired fermions. The effective ferromagnetic spin-spin coupling 

= 1, 2 N 



constants are given by = d^p^ for 



211 ] . In this sense, we may 



simply view the non-neutral charge bound state as a molecule with spin s, which could 
flip its spin antiparallel to the external fields H r due to thermal fluctuations. However, 
such spin fluctuations coupled to the channels of unpaired fermions and the spin-s charge 
bound states are suppressed by large fields H r at low temperatures. Thus the low energy 
physics is dominated by charge density fluctuations. Indeed we show that in the physically 
interesting regime T -C e r , A i+ u, 7 ^> 1 the breaking of charge bound states and spin wave 
fluctuations is strongly suppressed. Thus each dressed energy can be written in a single 
particle form e r (k) = h 2 rk 2 /2m — pS r > + 0(l/7 3 ), where the marginal scattering energies 
among composites and unpaired fermions and spin-wave thermal fluctuations are considered 
in the chemical potentials py>. 

The thermodynamics at finite temperatures can thus be obtained from the set of algebraic 
equations 
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r N Ap(i) 

u (r) w (r) _ y- y- tfL + fir) 

* * i{i + r-2j)\c\ +h ' 

i = j 



for r = 1, . . . , N. Here Li s (x) is the standard polylogarithm function and / s denotes the 
free energy of an effective ferromagnetic spin-spin interaction with coupling constant J^ r ) 



in the channel of the spin-s charge bound state [21] . However, because it is a spin singlet, 
there is no such effective spin-spin interaction for the spin-neutral bound state of iV-atoms, 
i.e., fs = and = 0. p^ r ) is the pressure for charge bound states of r-atoms. The 
pressure for the system is given by p = Y2^=i ■ The spin string contributions to thermal 
fluctuations in the physically interesting regime can be asymptotically calculated from the 
spin string equations for these non-spin-neutral charge bound states and unpaired fermions, 
i.e., fs^ ~ Te~ t t 1t where Io(x) = Y^T=o JWf(f ) 2k - ^ ^ s important to observe 

that /i r) becomes exponentially small as T — > 0. 

The suppression of spin fluctuations leads to a universality class of a multi-component 
TLL in each gapless phase, where the charge bound states of r-atoms form hard-core com- 
posite particles. In order to see this universal TLL physics, we calculate the leading low 
temperature corrections to the free energy by using Sommerfeld expansion with the TBA 
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equations to give 



/w/o- 



ttT 2 
~6?T 




(11) 



This result is consistent with the finite-size correction result (jl]). In the above equation 
/o = -Eg — J2^=i n rH r . This result proves the existence of TLL phases in ID gapped 
systems at low temperatures. Although there is no quantum phase transition in ID many- 
body systems at finite temperatures due to thermal fluctuations, the existence of the TLL 
leads to a crossover from relativistic dispersion to nonrelativistic dispersion between different 
regimes at low temperatures [19]. Nevertheless, we find that such a field-induced multi- 
component TLL only lies in a small portion of Zeeman parameter space. Linear Zeeman 
splitting may result in a two-component Luttinger liquid in a large portion of Zeeman 
parameter space at low temperatures. 

The thermodynamics of the gapless phase for the model ([1]) can be analytically calculated 
with linear Zeeman splitting, i.e. with A a+la = H for a = 1, . . . , N — 1. In the regime 
H c i < H < H c2 the TBA equations ( TTOT) reduce to two coupled equations for pS x > and p,( N >. 
The rest of the effective p,v) = 0. From these two equations, we can obtain the density 
n = dp/dfi, magnetization m z = dp/ OH and free energy / = fin — p by iteration. Fig. [2] 
shows the magnetization for attractive three-component fermions with a pure Zeeman field 
at different temperatures. We see that the linear field-dependent phase transitions in the 
vicinity of the critical points are smeared. 

The analytic results ( fTOl) and (lllj) indicate that the magnetization develops minima at 
the same temperatures as when the two-component TLL is broken. Indeed, a deviation from 
the linear temperature-dependent specific heat c v = ^ (-^ + and the magnetization 
minima occur around the same crossover temperatures - see the white-diamond-lines in Fig. 
[3] where we show a contour plot of the specific heat for SU(3) and SU(<i). The hard-core 
iV-atom composite particles and single atoms form an asymmetric two-component TLL lying 
below the white-diamond-lines. The filled-black-circles which separate the different regimes 
are determined by the magnetization values m z = and m z = 1. The peaks at the critical 
points are expected when the unpaired (bound state) band starts to fill (empty) [221 ]. 

To conclude, our analytic results for universal thermodynamics and quantum phase tran- 
sitions provide a unified description of attractive ID SU(N) fermions in the presence of 
external fields. Our formulae also provide the essential equation of states for studying 
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trapped A^-component Fermi gases with population imbalance. They thus pave the way for 
the study of multi-component TLL phases in ID systems of ultracold fermionic atoms. 
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(A) (B) 




FIG. 1: Phase transitions from the bound states of iV-fermions into (A) normal Fermi liquid and 
(B) paired states are induced by linear and nonlinear Zeeman splittings, respectively. Ellipses 
denote the charge bound states. 
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FIG. 2: Typical magnetization vs external field for SU(3) fermions with strong coupling (c = —10 
and n = 1) at different temperatures. The curves for SU(N) are similar. 
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(a)N=3,n=1.5 (b)N=4,n=2 




FIG. 3: Universal low temperature phase diagram: specific heat in the T-H plane for SU(3) and 
SU(4). H is rescaled by H° = 2cn/(N — 1). An asymmetric two-component TLL remains within 
the regime between H c \ < H < H C 2 below the white-diamond-lines. TLL of spin-neutral bound 
states of iV-atoms and TTL of fully-polarized fermionic atoms lie below the left and right filled 
black-circles, respectively. 
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